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ABSTRACT 


A three-dimensional, wire matrix, capacitance density meter was used to measure 
the density of liquid and two-phase hydrogen flow. The measurements were made in a 
4-inch (10. 16-cm) nominal diameter propellant feedline of a full-scale, simulated, 
nuclear -rocket test system. An evaluation of this capacitance density meter was made 
by comparing the density meter measurements with corresponding calculated densities 
in the feedline. The local hydrogen densities measured with the density meter were in 
relatively close agreement with the calculated densities. For the most part, the calcu- 
lated and measured densities differed by less than ±15 percent of the full-scale (liquid 
hydrogen) density. 
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SUMMARY 

A three-dimensional, wire matrix, capacitance density meter was used to measure 
the density of two-phase and liquid-hydrogen flow. The measurements were made in a 
4-inch (10.16-cm) nominal diameter propellant feedline of a full-scale, simulated, 
nuclear-rocket test system. An evaluation of this capacitance density meter was made 
by comparing the density meter measurements with corresponding calculated densities in 
the feedline. The results of two typical test runs are presented in this report. 

For the most part, the local hydrogen densities measured with the density meter 
were in relatively close agreement with the corresponding calculated densities. Most of 
the calculated and measured density values differed by less than ±15 percent of the full- 
scale (liquid hydrogen) density. It is concluded that this type of density meter could be 
used to measure two-phase and liquid-hydrogen densities with an accuracy that would be 
acceptable in many engineering applications. 


INTRODUCTION 

The measurement of fluid density or vapor content in flowing cryogenic systems is a 
subject of widespred interest. In the all-gas or all-liquid state, fluid density may be de- 
termined by measuring the fluid temperature and pressure and applying these measure- 
ments to an equation of state. This technique, however, cannot be used in a two-phase 
fluid region. 

In systems that contain two -phase flow, such as the nuclear rocket, the density of 
the flow mixture has a significant influence on heat-transfer rates and pressure drops in 
the system. And, in order to analyze such two-phase flow systems, it is necessary that 
the local fluid densities be known. 



Interest in the problem of two-phase density determination has led to the design and 
development of various instruments for measurement of density in croygenic flow sys- 
tems. The more promising of the many techniques proposed are based on either 
(1) measurement of the nuclear radiation attenuation properties of the fluid or (2) meas- 
urement of the average dielectric constant or capacitance of the fluid. In theory, both 
these measurable quantities are related to the fluid density. 

At the Lewis Research Center, considerable experience in density measurement has 
been gained from tests in which a capacitance density meter was used. In this report, 
results obtained from tests made with a meter using a three-dimensional, wire matrix, 
capacitance transducer are presented. A detailed description of this density meter is 
given in the section, DISCUSSION OF CAPACITANCE DENSITY METER. 

The capacitance transducer of this density meter was installed in the downstream 
end of the propellant transfer line (feedline) of a full-scale simulated nuclear-rocket test 
system. The meter was used to measure the local density of hydrogen in the propellant 
feedline during cold-flow startup experiments. The simulated nuclear-rocket test sys- 
tem in which this meter was used is located at the Plum Brook Station of the Lewis 
Research Center. (Ref. 1 has a detailed description of the system. ) 

The propellant feedline of the test system was instrumented so that the rate of heat 
transfer to the hydrogen flow during startup could be calculated directly from experi- 
mental measurements. The local hydrogen densities were computed by numerically in- 
tegrating the continuity and energy equations along the length of the propellant feedline. 
The calculated heat inputs, measured pressures, and measured inlet flow rates of the 
feedline were used in the density calculations. 

An evaluation of the capacitance density meter was made by comparing the density 
meter measurements with calculated local densities in the feedline. The measured and 
calculated local densities for two startup tests (runs 8 and 15) are presented and com- 
pared in this report. 


SYMBOLS 

2 2 

A free flow area, ft ; cm 

& coefficient in dielectric coefficient equation (eq. (5)), pF 

Jh rp dielectric coefficient (temperature dependent) (eq. (4)), pF 

c _ 

B capacitance temperature coefficient due to expansion of matrix (eq. (5)), pF/°R; 

pF/K 

C total capacitance, pF 
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Cp specific heat of fluid at constant pressure, Btu/(lb i T) & s g ) (°R) ; j/(g)(K) 

C g stray capacitance, pF 

C re f total capacitance at a particular reference condition, pF 
C e capacitance due to dielectric fluid in matrix, pF 

Cq capacitance of evacuated matrix assembly (eq. (8)), pF 

c heat capacity of stainless steel feedline, Btu/(lb mo __)( 0 R); J/(g)(K) 

D (C - c ref ) - <c 0 - C ref ) 

h enthalpy of fluid, Btu/lb masg ; j/g 

J mechanical equivalent of heat equal to 778 (ft)(lb force )/Btu (100 N-cm/J) 

K conversion factor (eq. (15)), equal to 144 in. 2 /ft 2 (10~^ mVcm 2 ) 

L length of flow element (eq. (A 14)) ft; cm 

M mass of feedline (eq. (23)), lb„_ 0 : g 

mass 7 

^ molecular weight, g/(g)(mole) 

m fluid mass, lt> mass ; g 

23 

N Avogadro’s number, 6. 0£<10 molecules/(g)(mole) 

p p 

P static pressure, lkf orce /in. abs; N/m 

3 3 

p specific polarization, ft /lb__ Q _ a ; cm /g 

mass 

Q rate of heat input to fluid per unit length of flow passage, Btu/(sec)(ft); 

j/(cm)(sec) 

Qrp rate of heat input to fluid from feedline inlet to feedline outlet (eq. (A20)) 

Btu/sec; J/sec 

T temperature, °R; K 

T temperature of matrix, °R; K 

V volume of a flow element (eq. (A 13)) ft 2 ; cm 2 

v velocity, ft/sec; cm/sec 

X ratio of mass flow rate of gas phase to total mass flow rate (eq. (A10)) 

x distance or position along flow path, ft; cm 

a molecular effective moment induced per unit field strength, cm /molecules 

e fraction of total liquid phase that is not converted to gas because of nonequilib- 

rium (eq. (All)) 
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e 


p 

r 

T 0 

d> 

Subscripts 

g 


j 

l 

M 

noneq 

o 

P o 

ref 

rms 

s 

sg 

si 

superheated 

gas 

w 

x 

(x-Ax) 

P 1 


dielectric constant 
density, lb mags /ft 3 ; g/cm 3 
time, sec 
zero time 

flow rate, lt> mass /sec; g/sec 
gas phase 

enthalpy of fluid at feedline inlet 

enthalpy of fluid at feedline outlet (feedline outlet corresponds to capaci- 
tance transducer inlet) 

conditions at feedline inlet 

j** 1 segment of flow passage (eq. (23)) 

liquid phase 

local average temperature of feedline 
nonequilibrium 

conditions at feedline outlet (or capacitance transducer inlet) 

static pressure at feedline outlet (feedline outlet corresponds to capaci- 
tance transducer inlet) 

reference condition (usually at ambient temperature and a pressure of 
approximately 1 psia (6. 9X1 0 3 N/m 3 )) 

root mean square 

static or stop action point of view (eq. (A 13)) 
saturated gas 
saturated liquid 

gas at temperature above saturation temperature (eq. (A9)) 

outside surface (wall) of feedline 

position x 

position (x-Ax) 

fluid density at condition 1 
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P 2 

fluid density at condition 2 

2 cp 

two-phase 

T 

time r 

(t-At) 

time (t-At) 

1 

condition 1 

2 

condition 2 


DENSITY METER 

The density meter described herein is composed of two major components: a three 
dimensional stainless -steel wire matrix capacitance transducer, and a capacitance-to- 
direct-current voltage converter. These two major components are described in the fol- 
lowing subsections. 



Figure 1. - Capacitance transducer. 


Electrical co nnectors 
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Description of Capacitance Transducer 


A photograph of the capacitance transducer is shown in figure 1. This wire matrix 
capacitor was designed and built by Space Science, Inc. (now Ikor Inc. , Burlington, 
Mass.). A sketch of the wire matrix capacitor is shown in figure 2. 

The assembly (fig. 2) consists of two electrically isolated three-dimensional wire 
matrices, each centered in the lattice of the other. (Solid and dashed lines are used in 
fig. 2 to represent the respective matrices.) The spacing between lattices is held fixed 
by a Teflon retainer ring. The three-dimensional matrix is mounted inside a pipe section 
and the electrical leads are brought through fittings on the pipe wall. 



Figure 2. - Matrix assembly of capacitance - density meter. (All dimensions are in inches (cm).) 
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The matrix capacitor senses the average dielectric constant of the fluid flowing 
through the matrix volume. The matrix assembly shown in figure 2 has approximately 
3200 cubical subvolumes, each with a side length of about 0. 1625 inch (0. 413 cm). Each 
of the subvolumes is equally sensitive to the fluid dielectric constant of its included 
mass. This is true regardless of the relative position of the subvolume in the overall 
sensed volume. This characteristic makes the matrix sensor insensitive to fluid inhomo- 
geneities, bubble position, or zero gravity. 

The density measured by the density meter is a static or M stop-action" density; that 
is, the measured density is the ratio of mass of fluid contained in the matrix to the vol- 
ume of the matrix. In the appendix, it is shown that the static or stop-action density of 
a two -phase flow mixture is equivalent to the so-called "mixing- cup" density when the 
velocities of the gas and liquid phases are equal. The mixing-cup density is defined as 
the ratio of total mass flow rate to the total volumetric flow rate. 


Description of Electronics 

The electronics for this density meter consists of a standard type of capacitance to 
direct-current (dc) voltage converter. A simplified illustration of a typical density 
meter installation is shown in figure 3. The dc voltage converter contains a radio- 
frequency (RF) oscillator, which is a transformer coupled to an inductance-capacitance 
(LC) bridge. 

The inductances of the bridge consist of autotransformers for coarse zero and span 
adjustments. One of the two capacitance legs includes the matrix capacitor. The other 



Figure 3. - Simplified illustration of density meter installation. 
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capacitance leg is a reference capacitor which is adjusted to balance the matrix capaci- 
tance. The output of the bridge is a signal proportional to the unbalance between the ma- 
trix and the reference capacitances. 

The capacitance to dc voltage converter was modified somewhat to accept a remotely 
acutated electrical calibration signal by substituting a predetermined electrical calibra- 
tion capacitor in place of the matrix capacitance. 

The density meter tested at Lewis had the following characteristics: 

O 

(1) Range, 0 to 5 pounds mass per cubic foot (0 to 0. 008 g/cm ) in liquid hydrogen 

(2) Accuracy, ±0. 10 pound mass per cubic foot (±0. 00016 g/cm'*) in liquid hydrogen 

(3) Hydrodynamic response time, 0. 63 (volume matrix) p/d) 

(4) Electronic response time, <0.01 second 

(5) Temperature effect, 2-percent shift in capacitance from 540° to 36° R (300° 

to 20° K) 

(6) Zero capacitance, 130 picofarads 

(7) Delta capacitance, 30 picofarads (gas to liquid hydrogen) 

(8) Transducer pipe size, 4 inch (10 cm) nominal schedule 40 

(9) Flow blockage, 13 percent of transducer pipe area 


Density Equation 

An equation relating density to the dielectric constant of a substance can be derived 
from classical electrostatic theory (see refs. 2 and 3). The equation, commonly re- 
ferred to as the Clausius -Mossotti equation, is as follows. 


NcA £ + 2 

.3 uT / 


( 1 ) 


The quantity inside the parentheses in equation (1) is called the specific polarization. 
For a particular substance, if oi is constant, the specific polarization of the substance 
is constant, and the substance is said to obey the Clausius-Mossotti equation. Equa- 
tion (1) was derived making use of the following assumptions: 

(1) Molecules are nonpolar and spherical. 

(2) No short-range interactions occur between molecules. 

(3) Molecules have isotropic spatial distribution. 

Although hydrogen meets none of these assumptions exactly over the entire range of 
densities, investigators have found that the specific polarization of hydrogen is essen- 
tially constant. Table I lists the values of specific polarization for hydrogen obtained by 
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TABLE I. - VALUES OF SPECIFIC POLARIZATION OF HYDROGEN 


Researchers 

Reference 

Hydrogen 

Temperatures, 

Approximate range 

Specific 



state 

K 

of specific 

polarization, 






volumes, 

% 






cm 3 /g 

cm /g 

Werner and Keesom 

8 

Liquid 

14 to 20.4 

13 to 14 

1.010*0.002 

Guillien 

9 



14 to 20.4 

13 to 14 

. 9984*0. 002 

van Itterbeek and Spaepen 

10 



20.4 

14 

. 988±0. 003 

Johns and Wilhelm 

11 



20.4 

14 

. 9986*0. 001 

Maryott and Buckley 

12 

Gas 

293 

12 000 

1.0104*0. 002 

Michels, et al. 

13 

Gas 

298 and 373 

16 to 11 000 

1.001*0.002 


different investigators. Regardless of the state of the hydrogen, the specific polarization 
(table I) equalled 1. 00 cubic centimeter per gram to within ±1 percent. 

Hust, et al. , (ref. 4) show that the dielectric constant of hydrogen depends only on 
fluid density, regardless of fluid phase. They report the dielectric constant of para- 
hydrogen for the range 24 to 100 K with densities from 0. 005 to 0. 080 gram per cubic 
centimeter. 

Therefore, the Clausius -Mossotti equation, which is sufficient to cover all three 
states of hydrogen, is 


p £ + 2 


( 2 ) 


where 


= constant =1.00 cmV g 

J / 

Equation (2) indicates that, for a constant specific polarization p, the fluid density is a 
unique function of the dielectric constant. And the average dielectric constant of the 
fluid is directly related to the capacitance. As a result of these relations, fluid density 
can be determined from capacitance measurements. 

The total capacitance of the three-dimensional assembly C is the sum of the stray 
capacitance C g and the capacitance of the dielectric fluid in the matrix Cg. Expressed 
as an equation, 



c = c s + c e 


( 3 ) 
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The capacitance due to the dielectric fluid in the matric is directly proportional 
to the dielectric constant of the fluid; that is, 

C £ = / T £ (4) 

The dielectric coefficient is a function of the matrix dimensions. And these 

A c 

dimensions change slightly with the matrix temperature. For a constant coefficient of 

thermal expansion, the dielectric coefficient Jhm is linearly proportional to the tem- 

A c 

perature of the matrix. 


= (j/+ BT C ) (5) 

By substituting the right. sides of equations (4) and (5) into equation (3), the following ex- 
pression for total capacitance of the matrix assembly is obtained: 

C = (s/ + BT c )£ + C g (6) 

The coefficient and the temperature coefficient of the matrix B are constants 
that are determined from calibration tests. The methods of determining these coeffi- 
cients are described later in this section. 

Capacitance is a difficult quantity to measure. Our experience has been that stray 
capacitance is always present in undeterminable amounts during calibrations, checkouts, 
and tests and that it always contributes an uncertainty to any measurement that is made. 
For the experiments descripted herein, the stray capacitance varied as much as ±4 pico- 
farads between tests. For hydrogen, the change in capacitance from the all gas to all 
liquid state is only about 30 picofarads. Consequently, the stray capacitance effect is 
significant. 

In order to overcome this problem, a technique was derived whereby differences in 
capacitances were measured instead of absolute values. A reference condition must be 
known during the data acquistion period. Usually, this is a known single-phase point or, 
ideally, a point at vacuum conditions. The differences in capacitance from this refer- 
ence capacitance are then measured. This technique successfully reduced stray capaci- 
tance measurement errors. To describe this technique, equation (2) is rewritten in the 
following equivalent form: 


p = I <E-D 
p (£ - 1) + 3 


( 7 ) 
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At perfect vacuum conditions, the dielectric constant £ is one. And the capacitance 
of the assembly at perfect vacuum conditions Cq is given by 


C 0 = (rf+ BT C ) + C s 


( 8 ) 


If the capacitances C and Cq are measured at a fixed matrix temperature T c 
equations (6) and (8) can be combined to give the following: 


£ = 


C - C 


s 



0 ) 


Subtracting one from each side of equation (9) gives 


£ - 


C - C 


1 = 


0 


C 0- C s 


( 10 ) 


Substituting this result into equation (7) and simplifying results in 

p = 1 < c - C Q> 

p (C - Cq) + 3(^+ BT C ) 


(U) 


Equation (11) may be used to evaluate density if both C and Cq are measured at the 
same temperature, T . 

For the tests reported herein, it was impossible to obtain perfect vacuum conditions 
in the experimental setup. As a result, the absolute value of the vacuum capacitance 
Cq could not be determined. 

In order to obtain an equation that could be used to evaluate density, equation (11) 

was expressed in terms of capacitance differences, which could be determined. A known 

single-phase reference condition existed during the data acquisition period just before 

the test run. During this period, the matrix was near ambient temperature and the 

3 2 

pressure was approximately 1 psia (6. 9x10 N/m ). The measured capacitance at this 
reference condition of low pressure and ambient temperature C re ^ can be expressed in 
the following form: 


C ref “ B1 W £ ref + C s 


( 12 ) 
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Substracting C rg £ from both C and Cq in equation (11) and substituting 

for j/+ BT„ results in an equation for density in terms of differences in capacitances: 
c 


p 


1 

p 


(c 


C ref> - < C 0 ‘ C ref> 


(c - C ref ) - (Cq - C ref ) + 3 ^ t 


(13) 


As will be shown, the capacitance difference terms and the coefficient in equa- 

**■ p 

tion (13) can be readily evaluated. 

The quantity C - C re j is the change in capacitance compared with the low pressure, 
ambient temperature reference condition existing before a test run. This is the quantity 
that is experimentally measured. The quantity Cq - C rg j is a correction term to ac- 
count for (1) the change in the matrix dimensions due to temperature changes and (2) the 
nonvacuum zero conditions. This correction term is given by the following equation: 


(Cq - C ref ) = 1 - £ ref ) + £ ref B(T c - T ref ) 


(14) 


The coefficients A— and B in equation (14) are given by 
A c 



(15) 


B - f C 0\ J%1 ' C 0,2 


v3T 


'p- 0 


T l- T 2 


f p = 0 


(16) 


The density measurement system must be calibrated to obtain numerical values for 
Jhrr, and B. To determine the dielectric coefficient Jj - ) T , the capacitance is measured 

at two different known densities (p^ and Pg) with the matrix at constant temperature T c . 

The dielectric constants (£ and £ ) in equation (15) are computed by substituting the 

\ Pi Ptf 

known densities p^ and Pg into the Clausius-Mossotti equation. The numerical values 

of C C , £ , and £ are then substituted into equation (15) to evaluate /!„ . 

P 1 @2 P 1 p 2 A c 

To determine B, the capacitance of the matrix assembly is measured at two known 

temperatures (Tj and Tg) with the matrix in a vacuum. The measured capacitances, 

Cq i and Cq ^ along with the known temperatures are then used in equation (16) to 
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evaluate the coefficient B. The density of the fluid is then obtained by substituting the 

numerical values of Aj, , ( C - C ref ), and (C Q - C ref ) into equation (13). By an inde- 

c 

pendent calibration, the value of was determined to be 131. 82 picofarads and the 

c 

value of B was 0. 00235 picofarad pel* R (0. 00423 pF/K). 

In the appendix, the accuracy of this particular density meter is examined. And an 
rms inherent error for this meter is calculated to be about 3 percent (full scale) when 
used to measure hydrogen densities. This 3 -percent error is based on the estimated un- 
certainties in each of the terms in equation (13). 


RESEARCH APPARATUS AND PROCEDURE 

The density meter was tested in a full-scale simulated nuclear-rocket test system. 
The full-scale test system consisted of a liquid-hydrogen storage tank, turbopump, pro- 
pellant feedline, reactor, and nozzle assembly. The assembled configuration of this test 
system is shown in figure 4. Reference 1 provides a detailed description of the major 
components of this test system. 

In this report, the section of propellant feedline which extends from the main dis- 
charge valve to the capacitance transducer is of major importance. The relative posi- 
tions of the main discharge valve and the capacitance transducer are shown in figure 4. 


Propellant Feedline 

The feedline considered in this study is a 4-inch (10-cm) nominal diameter, sched- 
ule 5, AISI 304, stainless-steel pipe. This pipe has an inside diameter of 4. 334 inches 
(11. 0 cm), a wall thickness of 0. 084 inch (0. 213 cm), and an overall length (from main 
valve to capacitance density meter) of approximately 17 feet (5.2 m). 

A sketch of the feedline is shown in figure 5. The capacitance transducer is shown 
near the discharge end of the feedline in figure 5. Two gimbal expansion joints were in- 
stalled in the vertical leg of the feedline to allow movement of the line during chilldown. 
The entire feedline was covered with a 4-inch (10-cm) thick layer of polyurethane foam 
insulation. The purpose of the insulation was to minimize heat transfer to the line by 
free convection and to eliminate possible condensation of air on the outside surface of 
the pipe. 
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Figure 4. - Full-scale, cold-flow nuclear-rocket test system. 



s^-Tu rbine flowmeter 

\ 



Figures. - Propellant feedline instrumentation. 




Feedline Instrumentation and Accuracy 


The locations of the pressure and material temperature measurements for the feed- 
line are shown in figure 5. Pressures were measured at four points on the feedline with 
calibrated strain-gage transducers. The full-scale range of the transducers was 
200 psia (13. 8x10® N/m^). Temperatures of the stainless -steel feedline were measured 
with fine wire (30 gage or 0. 0255 cm diam) copper -constantan thermocouples. The 
measuring junctions of these thermocouples were welded to the outside wall of the feed- 
line. 

The hydrogen flow rate at the feedline inlet was measured with a turbine flowmeter. 
The turbine meter is shown in figure 5 just upstream of the main discharge valve. 

Besides the material temperature and pressure-measurement stations listed in fig- 
ure 5, other stations for the measurement of fluid temperatures and pressures were just 
ahead of the main discharge valve. These fluid temperatures and pressures were used 
to determine the enthalpy and density of the liquid hydrogen upstream of the main dis- 
charge valve. The fluid temperatures were measured with precision resistance sensors. 

The estimated accuracies of the different experimental measurements are listed in 
table II. The accuracies listed for copper-constantan thermocouples in the range from 
200° to -300° F (660° to 160° R or 367 to 89 K) were taken from the National Bureau of 
Standards specifications (ref. 5). For temperatures below -300° F (160° R or 89 K) the 
accuracy of measurements made with copper-constantan thermocouples is difficult to 
access. In the range from -300° to -430° F (160° to 30° R or 89 to 16. 7 K), the sensi- 
tivity of copper-constantan is small. (Sensitivity as used herein is the change in thermo- 
electric potential with respect to temperature.) Because of the small sensitivity, meas- 
urements made with copper-constantan near the liquid-hydrogen temperature regime are 
appreciably affected by factors such as (1) resolution of the data recording system, 


TABLE n. - INSTRUMENTATION ACCURACY 


Variable 

Instrument 

Estimated accuracy 

Flow rate 

Turbine flowmeter 

±2 percent of actual flow rate 

Pressure 

Strain- gage transducer (0 to 200 psia) 

±1/2 percent of full scale, i. e. , ±1 psi 


(0 to 13. 8x1 0 5 N/m 2 ) 

(±6. 9X10 3 N/m 2 ) 

Fluid temperature 

Resistance-type sensor 

±0. 50° R (±0. 28 K) in temperature 
range from 30° to 60° R 
(17 - 33 K) 

Material temperature 

Copper-constantan thermocouple 

±0. 75° F (±0. 42 K) for temperatures 
ranging from 200° to -75° F (367 to 
to 214 K); ±1 percent of Fahrenheit 
temperature from -75° to -300° F 
(214 to 89 K) 
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(2) variations in the thermoelectric powers of different wire pairs, and (3) line noise in the 
data acquisition system. For the tests considered in this report, the errors associated 
with material temperature measurements made in the liquid-hydrogen temperature re- 
gion were estimated to be as much as ±10° R (±5. 6 K). 

Test Operations 

The experimental procedure used in the operation of the full-scale, simulated 
nuclear-rocket test system is described in reference 1. A brief summary of the proc- 
esses occurring in the propellant feedline during a typical startup test is given below. 

For the test runs described herein, the turbopump and piping upstream of the main 
discharge valve were precooled and filled with liquid hydrogen before the start of a test. 
The temperatures of the turbopump and piping upstream of the main valve were reduced 
to approximately the temperature of liquid hydrogen. The temperatures of the feedline 
downstream of the main discharge valve were initially near ambient temperature. 

The flow of hydrogen in the feedline is initiated by a programmed opening of the 
main discharge valve. Time zero Tq corresponds to the initiation of flow by opening 
the main discharge valve. 

In the early part of a startup test run, the feedline is near ambient temperature, 
and a large fraction of the liquid-hydrogen flow entering the feedline is changed to vapor 
as it moves downstream in the line. The hydrogen flow arriving at the capacitance 
transducer is therefore predominately in the vapor state. As time increases, the turbo- 
pump causes an increase in the flow and pressure in the feedline; concurrently, the tem- 
peratures of the propellant feedline decrease. As a result, less heat is transferred 
from the feedline to the fluid, and the density of the hydrogen arriving at the capacitance 
transducer increases as the rim continues. Eventually, the flow exists as a subcooled 
liquid phase throughout the feedline. 

The density meter was used to measure liquid and two-phase hydrogen densities in 
more than 30 bootstrap startup experiments. The range of variables covered in all the 
bootstrap tests, that is, the pressures, fluid temperatures, and flow rates, were nearly 
the same. The major difference between these bootstrap tests was in the flow- rate -time 
profiles. In some test runs, the feedline inlet flow rates were relatively smooth; in 
other runs, however, they were oscillatory, especially near the beginning of the runs. 
The two test runs presented in this report represent one of each of the two classes of 
flow-rate - time profiles. Run 15 is typical of those runs that had a smooth flow-rate 
profile, and rim 8 is typical of those in which the flow-rate - time profile was oscillatory. 

The two runs (8 and 15) were analyzed in considerable detail. Most of the other test 
runs were also analyzed. However, the calculational procedure used was somewhat sim- 
plified. The calculated densities obtained for these other tests were used primarily as 
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TABLE ID. - TEST VARIABLES 


Run 

Length 

Flow rate range 

Pressure range at feedline inlet 

Fluid temperature at feedline inlet 


of test, 
sec 

lb /sec 

mass' 

kg/sec 

psia 

kN/m 2 

°R 

K 

8 

18 

0 to 32. 5 

0 to 14. 8 

0. 5 to 172 

3.45 to 1.19X10 3 

39.4 to 45.3 

21.9 to 25.2 

15 

16 

0 to 31.5 

0 to 14. 3 

0. 5 to 160 

3.45 to 1.10X10 3 

39.2 to 42.0 

21.8 to 23.3 


qualitative information. This information was considered, however, in the evaluation of 
the density meter. 

For both tests (runs 8 and 15) the liquid-hydrogen storage tank was pressurized to 
35 psia (2. 42xl0 5 N/m 2 ). The pressures in the feedline downstream of the main valve 
at Tq were approximately 1. 0 psia (6. 9x10 N/m ). The important variables for runs 
8 and 15 are listed in table m. The feedline inlet conditions for the two runs are pre- 
sented in the section RESULTS AND DISCUSSION. 

Data Acquisition and Processing 

All temperature and pressure measurements in the propellant feedline were re- 
corded on a digital multiplex system. The recorded data were then processed on a digi- 
tal computer, which converted the recorded millivolt values to engineering units of tem- 
perature and pressure. A description of the data acquistion system is given in refer- 
ence 1. 


ANALYTICAL PROCEDURE 

Consider the case of a one- or two-phase homogeneous fluid flowing in a constant 
area passage with heat addition. The one -dimensional differential equations of conti- 
nuity and energy for nonsteady flow are written as follows: 

Continuity: 


A ie + i“ = 0 

3r 0 x 

Energy: (neglecting kinetic energy term) 

Ap — - — + <i> — = Q* 

0T J 0T 0X 


(17) 


(18) 
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Now, consider a finite segment of passage with length Ax, such as that shown in 
figure 6. In figure 6 there are two sets of inlet and outlet conditions: one set corre- 
sponds to time t - At and the other set to time t. If the variable designations in fig- 
ure 6 are used and if a finite time interval At is considered, the differential equations 
of continuity and energy (eq. (17) and (18)) can be expressed in terms of finite differences. 
The corresponding equations are 

Continuity: 


A o Aw 
_J_+ _*_ = 0 
At Ax 


(19) 


Energy: 

A h AK A P A h 

A p — — + u) = Q (20) 

At J At Ax 

Expanding equations (19) and (20) and using appropriate inlet and outlet values to repre- 
sent averages, results in 

Continuity (expanded finite difference form): 

a R P x, t + P x-Ax, P ~ fos, t-At + P x-Ax, T-Ap j + ^x, t ~ ^x-Ax, t _ q ^21) 

L 2 At J Ax 
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Energy (expanded finite difference form): 


A <P X , r + P x-A x, ■? 

2 


K T + h x-Ax, P ~ fex, t-At + h x-Ax, t-At ^ 
2 At 


AK 

J 


^ P x, T + P x-Ax , 


/p + p V 

' x, r-Ar x-Ax, t-At’ 

2 At 


+ 


CO + CO A 

X, T x-A x, T 


JElZ. 


~ h x- Ax, 
Ax 
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The expanded finite difference forms of the continuity and energy equations were 
used to calculate the local fluid densities in the propellant feedline as a function of time. 
In applying these equations, the feedline (from main valve to capacitance transducer) was 
considered to be made up of 28 segments arranged in series. The average length of the 
segments was approximately 7. 5 inches (19. 05 cm). The inlet of the first segment cor- 
responds to the position x = 0 in figure 5. 

The heat input to each segment was determined from the stainless -steel feedline 
wall temperature measurements. As stated in the RESEARCH APPARATUS AND PRO- 
CEDURE section, wall temperatures were measured at 16 different stations along the 
length of the feedline. These temperature measurements were linearly interpolated, and 
an appropriate wall temperature was determined at each point in time for each of the 
28 segments. In the calculations, it was assumed that the measured local wall tempera- 
tures were equal to the average local material temperatures of the feedline; that is, the 
temperature gradient across the feedline wall thickness was assumed to be negligible. 
Hence, the heat input to the fluid in time interval At from wall segment j with length 
AXj is given by the following equation: 



A similar interpolation procedure was used to determine the pressures at each point 
in time for the endpoints of the 28 segments. 

Throughout the experimental runs, the hydrogen just ahead of the main discharge 
valve was subcooled (i. e. , the temperature of the liquid hydrogen was below the boiling 
point at the existing pressure). And the enthalpy there was determined by measuring the 
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local fluid pressure and temperature and applying these measurements to an equation of 
state. 

To determine the enthalpy at the inlet to the first segment, a constant enthalpy throt- 
tling process across the main discharge valve was assumed. Hence, at each point in 
time, the enthalpy at the inlet to the first segment was equal to the enthalpy just ahead 
of the main discharge valve. 

The hydrogen density at the inlet to the first segment was determined as a function 
of time from the local enthalpy and pressure values. The following procedure was used 
to get the local densities, enthalpies, and flow rates in the feedline during a flow test. 

Beginning with the first segment, at time zero, and considering the time interval 
At, the expanded finite difference forms of the continuity and energy equations were 
solved simultaneously. And the density, enthalpy, and flow rate at the outlet of the first 
segment at time t (i. e. , p , h and d> ) were obtained. These outlet conditions 

Xj T Xj T Xj T 

from the first segment were then assigned to the inlet of the following segment and the 
calculations were repeated. This stepwise procedure was continued and the fluid condi- 
tions at the outlet of each of the 28 segments at time t were determined. 

These calculations were repeated for each of the 28 segments and for each time in- 
terval At to the end of the flow test. From these calculations, the density of hydrogen 
at the inlet to the capacitance transducer was determined as a function of time. (The in- 
let of the capacitance transducer corresponds to the outlet of segment 28.) 

For the two test runs described in this report, a constant time interval of 0. 2 second 
was used. The calculations were performed on a high-speed digital computer. The hy- 
drogen properties used in the calculations were obtained from reference 6. Because of 
dependency of p on h an iterative procedure was required for the solution of 

Xj T Xj T 

the equations. 


RESULTS AND DISCUSSION 
Comparison of Measured and Calculated Densities 

Figure 7 show the measured and calculated hydrogen densities at the feedline outlet 
as a function of time for runs 8 and 15. The measured and calculated traces of density 
for each of these tests have a similar shape. And, for the most part, the hydrogen den- 
sities measured with the capacitance density meter are considered to be in good agree- 
ment with the calculated densities. 

In both of these tests, the hydrogen at the feedline outlet was a two-phase (gas-liquid) 
mixture for times from about 0. 5 to 4. 5 seconds. 
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0 2 4 6 8 1 0 12 14 16 18 


Time, t, sec 
(b) Run 15. 

Figure 7. - Measured and calculated hydrogen density at feedline outlet as func- 
tion of time. 

The accuracy of the measured and calculated densities are examined in the appen- 
dix. In addition to density meter error, four important items are considered that could 
have a bearing on the comparisons of measured and calculated two-phase densities: 

(1) Thermodynamic nonequilibrium between the gas and liquid phases 

(2) Velocity differences between gas and liquid phases 

(3) Possible inaccuracies in experimental measurements used in density calculations 

(4) Approximations used in the determination of heat-transfer rates 

The magnitude of the errors in two-phase density resulting from items (1) and (2) 
depends on the type of two-phase flow and the degree of turbulent mixing between the gas 
and liquid phases. For the test runs reported herein, no attempt was made to experi- 
mentally evaluate the temperature and velocity differences between the gas and liquid 
phases. Therefore, the affects of items (1) and (2) could not be quantatively accessed. 
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However, an indication of these individual effects on two-phase hydrogen density is given 
in the appendix. 

For a temperature difference of 10° R (5. 56 K) between the gas and liquid phases, 
the maximum relative error in calculated density is estimated to be about -20 percent. 

And for gas-to-liquid phase velocity ratios ranging from 0. 75 to 1. 25, the resulting 
maximum relative error in calculated two-phase density is approximately ±15 percent. 

The errors in calculated two-phase densities that may result from items (3) and (4) 
are of a systematic class. The magnitude of these possible errors are estimated in the 
appendix. 

The maximum relative error in calculated two-phase hydrogen density due to item 
(3) was estimated to be approximately ±9 percent. And the maximum relative error in 
calculated density due to item (4) is about ±5 percent. 

Most of the calculated and measured densities in figure 7 differed by less than 
±15 percent of the full-scale (liquid hydrogen) density. When the possible effects of the 
four items discussed previously are considered, the agreement in measured and calcu- 
lated densities (fig. 7) is good. For the most part the differences between calculated 
and measured two-phase densities are within the uncertainty band that results from the 
systematic errors, that is, items (3) and (4). 

An important feature of this density meter is its reliability. The capacitance trans- 
ducer proved to be a rugged instrument. It was installed in the propellant feedline and 
subjected to more than 30 transient cold-flow startup tests. After used in all these tests, 
there was no apparent damage to the transducer. 

In this report, data for two test runs in which the density meter was used are pre- 
sented. For the most part, the emphasis was placed on determining the relative accu- 
racy of the measured local densities for these tests. It was concluded that measure- 
ments made with this density meter were relatively accurate. Information obtained from 
qualitative comparisons of calculated and measured densities for other test runs support 
this conclusion. 


Input Data 

The time -dependent fluid conditions at the inlet to the propellant feedline for the two 
test runs are shown in figures 8 to 11. 

The hydrogen-inlet flow rates are shown as a function of time in figure 8. Figures 9 
and 10 show the inlet fluid temperatures and pressures against time. Fluid enthalpies at 
the feedline inlet are shown against time in figure 11. (The differences between the solid 
and dashed lines in fig. 11 are a measure of the subcooling during these test runs.) 


23 



I 



(b) Run 15. 

Figure 9. - Fluid temperature at feedline inlet against time. 


As discussed in the section ANALYTICAL PROCEDURE, these time varying inlet 
conditions were used in the finite difference equations of continuity and energy to calcu- 
late hydrogen densities at the feedline outlet. As stated in the RESEARCH APPARATUS 
AND PROCEDURE section, the major difference between the two reported runs was in 
the flow-rate - time profiles. In run 15 the flow-rate change with respect to time at the 
feedline inlet (fig. 8(b)) was relatively smooth. And the measured and calculated density 
traces for run 15 (fig. 7(b)) were also relatively smooth. In run 8, however, the feed- 
line inlet flow rate (fig. 8(a)) was oscillatory in the early part of the test; the corre- 
sponding density traces for rim 8 (fig. 7(a)) have a similar oscillatory behavior. 
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(b) Run 15. 

Figure 11. - Fluid enthalpy at feedline inlet against time. 



CONCLUSION 


For the rims presented herein, the local hydrogen densities measured with the ca- 
pacitance density meter were in relatively close agreement with the corresponding cal- 
culated densities. 

Most of the calculated and measured density values differed by less than ±15 percent 
of the full-scale (liquid hydrogen) density. For the most part, the calculated and meas- 
ured densities for these tests were within the bands of uncertainty that resulted from (1) 
the inaccuracy of measurements that were used in the analytical density calculations and 
(2) the inaccuracy of the density meter. 

It is concluded that the fluid density measurements made with this capacitance den- 
sity meter were relatively accurate. Furthermore, this type of density meter could be 
used to measure two-phase and liquid-hydrogen densities with an accuracy that would be 
acceptable in many engineering applications. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 11, 1968, 

122-29-03-04-22. 
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APPENDIX - ERROR CONSIDERATIONS 

The discrepancies between measured and calculated densities result from three dif- 
ferent sources: (1) the errors of the density meter itself, (2) the errors in the data and 
the assumptions used to determine the calculated density, and (3) deviation of the experi- 
ment from the assumed conditions of thermal equilibrium and zero slip between the fluid 
phases. Each of these error sources are analyzed in the following sections of this ap- 
pendix. 


Error of Density Meter 


The errors in the density meter correlation can be attributed to three different 
sources. The first is the ±1 percent uncertainty in the specific polarization in the 
Clausius -Mossotti equation. The second is the uncertainty in the calibration coefficient. 
The third and major source of error is the experimental uncertainty in measuring the 
matrix capacitance. 

When D = (C - C f ) - (C - C ref ) is substituted into the density equation (eq. (13)), 
the following error equation is derived: 


dp = 

dp 

1 

3iL 

L c dD 

L 

3 d jft 

P 

P 

T 

D + 3# t D 
c 


D + 3j^ t ^ T c 
c 


But 


3 3 

D+3 $ T 3J T 

c c 

is a good approximation because 3J ^ is much larger than D. Therefore, the root 
mean square (rms) error is l ' 
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The relative error in specific polarization p is about ±1 percent. Table I shows 
the consistency of p for six different experimenters. The ±l-percent variation in p 
shown in this table includes the random variation between experimenters and any nonran- 
dom variation in the specific polarization due to density dependency. 

The nonrandom variation is a result of p being a weak function of density. This 
variation indicates the degree to which hydrogen fluid does not obey the Clausius- 
Mossotti equation. Only the nonrandom varition in p will enter in the density meter 
equation error, because, here, p was assumed to be a constant. The individual calibra- 
tion of the meter (i. e. , the determination of the factor Jhrj, ) should compensate for the 

c 

experimenter’s random error. Table I indicates little if any nonrandom error. Most of 
the variation appears to be random variation among the experimenters. Hust, et al. , 

(ref. 10) report that the specific polarization is a function of density only. They give the 
relation between specific polarization (cm /gm) and density (gm/cm ) by 

1= 0.9957 - 0.09069 p+ 1.1227 p 2 (A2a) 

P 

From this relation, it can be shown that the maximum change in specific polarization is 
less than 0. 2 percent. This value will be used in the error analysis. 

The second error, the calibration uncertainty, appears as a variation in $ T , the 

dielectric coefficient of the matrix. This error comes from two main sources. The de- 
termination of the density of the calibration fluid has some uncertainty. Also, the meas- 
ure of the output capacitance of the system during calibration is affected by stray capaci- 
tance and instrumentation error. 

Starting with the calibration equation for the dielectric coefficient of the matrix 
(eq. (15)), the rms calibration error is 



But the error in the dielectric constant is related to the error in calibration density as 
follows: 


d£ = 2 (£ + 2) 2 dp 

3 


(A4) 
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where £ is an average value in the range of <£.. For simplicity, the approximation is 

£ = £„ 

p 2 

But 


dp = 



(A5) 



Reasonable estimates of calibration measurement errors are ±^ picofarad in capaci- 
tance, ±1 psia in pressure (±6. 9x10^ N/m^); and ±1° R (±5/9 K) in temperature. There- 
fore, the rms error due to the calibration is 2. 27 percent. 

The third error source is the uncertainty in measuring capacitance during an experi- 
mental run. This error is a result of changes in stray capacitance and also instrumen- 
tation error. The rms experimental error equation is 


a 


rms 


d(c - C ref ) 


((C - c ref ) - (C 



d(C f 


'ref' 


'ref 


•) - (C f 


'ref. 


1/2 


(A7) 


Estimates of experimental measurement errors are ±^ picofarad in d(Cg - C re j) and 
picofarad in d(C - C re j). Therefore, the rms error due to experimental measure- 
ment is 1. 87 percent. 

Total density measurement system rms error can now be found from equation (13). 
The total rms error is 2. 98 percent. And this is the rms error of the density meter 
only. 


Items That Affect Comparison of Calculated and Measured Two-Phase Densities 

In the following paragraphs, four important items that could affect the comparison of 
calculated and measured two-phase hydrogen densities are considered. The first and 
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second items discussed deal with the possible differences in the calculated and measured 
densities due to nonuniformities in the temperatures and velocities of the phases in a 
two-phase flow mixture. The third and fourth items deal with possible inaccuracies in 
the calculated densities resulting from errors in experimental measurements and ap- 
proximations used in the density calculations. 

Thermodynamic nonequ ilibrium between phases . - The calculated two-phase den- 
sities in figure 7 are based on the assumption that the gas and liquid phases at the feed- 
line outlet were in thermodynamic equilibrium. Heat-transfer experiments conducted 
with two-phase hydrogen flow(e. g. , ref. 7) have shown that temperature differences may 
exist between the gas and liquid phases. With nonequilibrium conditions, it can be shown 
that the corresponding average two-phase density is always less than the equilibrium 
two-phase density. 

For the test runs described in this report, no attempts were made to experimentally 
determine the degree of nonequilibrium in the two-phase flow. But it is likely that some 
differences did exist between the gas and liquid phase temperatures. 

The effect of nonequilibrium on two-phase density is illustrated in figure 12. The 
nonequilibrium values of density are based on an assumed 10° R (5. 56 K) superheat of 
the gas phase. 

The equilibrium two-phase densities in figure 12 were calculated from the following 
equation: 



Figure 12. - Equilibrium and nonequilibrium two-phase densities against quality for pres- 
sures of 35 psia (2.4M0 5 N/m 2 ) and 100 psia (6.9xl0 5 N/m 2 ). Nonequilibrium densities 
based on 10° R (5. 56 K) superheat of gas phase. 
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P 2cp~ x 


1 - X 


(A8) 


sg ^s l 


The nonequilibrium two-phase densities were computed from the equation 


2 (p, noneq 


X - e 


1 + e - X 


(A9) 


^superheated gas p sl 


The fluid quality X in equation (A8) and (A9) is defined as 


x _ Mass flow rate of gas phase _ ^g 


(A10) 


Total mass flow rate 


u) + Ct), 

g l 


In equation (A9), e is the fraction of the total liquid phase, which is not converted to 
gas because of nonequilibrium conditions. The value of e depends on the pressure 
level and the amount of superheat of the gas phase. At a fixed pressure, e is defined as 
follows: 


x f i 

*/at 


dT 


e = 


superheated gas 


(»Sg - V + 


f c p dT 

^ supei 


superheated gas 


(All) 


Figure 12 shows the equilibrium and nonequilibrium two-phase densities as a func- 
tion of quality X for pressures of 35 and 100 psia (2. 41X10 5 and 6. 9xl0 5 N/m 2 ). The 
nonequilibrium traces in figure 12 are based on a 10° R (5. 56 K) superheat of the gas 
phase. 

The relative difference between equilibrium and nonequilibrium densities, that is, 
(Pg ^ - P2cp noneq} / p 2 cp’ increases as the equilibrium two-phase density decreases. For 
the pressure range from 35 to 100 psia (2. 41X10 5 to 6. 9xl0 5 N/m 2 ) the maximum rela- 
tive difference between equilibrium and nonequilibrium densities due to a 10° R (5. 56 K) 
superheat of the gas phase is approximately -20 percent. And the maximum difference 
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expressed as a fraction of the full-scale (liquid hydrogen) density, that is, 

{p 2cp “ p 2cp, noneq) /p sZ’ is about " 6 P^cent. 

The difference term (pg^ - P^y none q) is approximately equal to the calculated den- 
sity minus the measured density. This difference results when the gas and liquid phases 
are at different temperatures. Therefore, with nonequilibrium two-phase flow, the cal- 
culated densities would be expected to be less than the corresponding measured densities. 

Velocity differences between phases . - The equilibrium density of a two-phase (gas- 
liquid) flow mixture is given by equations (A8) and (A 10). By combining these equations, 
the equilibrium two-phase density can also be expressed as 


Total mass flow rate _ ^g + w l 
^ Total volumetric flow rate co u, 

S- + -L 

p sg p sl 


(A 12) 


In this form, the density of a two-phase flow mixture is dependent only on the den- 
sities of the saturated gas and liquid phases and the flow rates of the phases. The den- 
sity given by equation (A 12) is commonly referred to as mixing- cup density. 

Now consider a finite element of constant area passage with two-phase flow as illus- 
trated in figure 13. From a static or stop action point of view, the average density of 
the flow mixture in the element of volume V (fig. 13) is given by 


P 


2<p,s 


Total mass in V 
V 


m g +m ; 

v 


(A 13) 


The mass of gas in volume element V is 


o;„L 


m„ = 


g 


(A 14) 


Flow 


Figure 13. - Element of constant area passage with two-phase flow. 
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And the mass of liquid in volume element V is 


m 


l 





(A 15) 


The volume of the element is 


m m d>„L 
V = — + — — = 




p sg p sl v g p sg v l p sl 


(A 16) 


Substituting the right side of equations (A 14), (A 15), and (A 16) into equation (A 13) and 
simplifying gives 


p 2cp,s 



(A 17) 


The density given by equation (A 17) is the theoretical density of the two-phase fluid 
as measured by the density meter. The ratio v g/ v j equation (A 17) is commonly re- 
ferred to as slip. For slip equal to one, equation (A17) is identical to equation (A12). 

For Vg /vj # 1, the density given by equation (A17) is different from that given by equa- 
tion (A 12), and kP2(p ~ p 2<p s^ is the difference between calculated and measured two- 
phase densities due to nonuniform phase velocities. 

In general, the average velocities of the individual phases of a two-phase flow mix- 
ture are not the same. For steady flow in a horizontal passage with heat addition, the 
slip Vg/v^ is usually greater than one. For the experiments described in this report, 
the flow was unsteady. And the slip at the capacitance transducer may have been greater 
than one at some times during these tests and less than one at other times. 

The affect of slip on the density of two-phase hydrogen is shown in figure 14. (The 
curves in fig. 14 are based on an assumed constant pressure of 40 psia (2. 76X10^ 

N/m ).) For slips ranging from 0. 75 to 1. 25, the maximum relative difference in meas- 

pr p 

ured and calculated densities at 40 psia (2. 76x10 N/m ), that is, i.P2cp ~ p 2cp s^ p 2cp’ 

is approximately ±15 percent. And the maximum difference as a fraction of full-scale 
(liquid hydrogen) density (P 2 ^ " ^2cp S )/P S i * s a ^ ou ^ ±5 percent. At pressures greater 
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Figure 14. - Two-phase hydrogen density against quality with slip as parameter. (Data 
based on pressure of 40 psia (2. 76xl0 5 N/m 2 ). ) 

than 40 psia (2. 76X10 5 N/m 2 ), the maximum difference is somewhat less than the previ- 
ously stated values. 

Accuracy of experimental measurements used in density calculations . - The esti- 
mated accuracies of the measured flow rates, pressures, fluid temperatures, and ma- 
terial temperatures are given in table n. Each of these measurements was used in the 
calculations to determine the density of hydrogen at the feedline outlet. In the discussion 
that follows, the effect of errors in these measurements on the calculated two-phase 
density at the feedline outlet is examined. 

In the two-phase region, the local fluid density can be expressed as a function of the 
local pressure and enthalpy; that is, 


'V =( V P o- h o > 


(A18) 


And the total differential of the two-phase density is given by 



(A 19) 


If the time derivatives in equation (20) are neglected and uniform, steady flow is as- 
sumed, the enthalpy at the feedline outlet h Q can be written as 
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(A20) 


J Cx 

L Q 


Q dx = hj + 


Qn 


CO 


The total differential of the fluid enthalpy at the feedline outlet is 



The partial derivations in (A21) can also be written as 



(A21) 


(A22) 


(A23) 


(A24) 


If the right sides of equations (A21) to (A24) are substituted into (A 19), the relative 
error in calculated two-phase density is given as 

r 



p 2 cp p 2cp 



o' U 


Q t deb 


(coy 


K 


<Jt) 


M 


(A25) 


The relative error given by (A25) is dependent on the numerical values of P Q , Pg^, Qj>, 
and c b. And for fixed values of P Q , Q T , and cb, the relative error in calculated two- 
phase density approaches a maximum as p 2( ^ approaches the saturated liquid density. 
For P = 50 psia (3. 45x10^ N/m^), co = 20 pounds mass per second (9. 1 kg/sec), 

• O c 

Qrj, = 100 Btu per second (1. 055x10 j/sec), and p 2 ^, = 4 pounds mass per cubic foot 
(0. 064 g/cm 3 ), the approximate numerical values of the partial derivations in equa- 
tion (A25) are 
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0. 025 


< lb mass><‘"- j 
( Ib force)< tt3 > 



(g)^ 2 ) 

(N)(cm 3 ) 



-0.32 


lb 


mass 


(Btu) (ft 



And from table II, 


dP Q = ±2 psi (1. 38X10 4 N/m 2 ) 


^ = ± 0 . 02 

u> 


From the specified accuracy of the fluid temperature measurements in table II, dh. is 
about ±0. 75 Btu per pound mass (±1. 76X10 J/kg). And, from the specified accuracy of 
the material temperature measurements, dQrp was estimated to be about 3 Btu per sec- 
ond (3. 165X10 3 j/sec). 

Using these numerical values in equation (A25) results in a maximum relative error 
in P 2 (p> due to possible inaccuracies in experimental measurements, of about 9 percent; 
that is, dp 2(p /p 2(p s 0. 09. 

Approximations used in determination of heat transfer rates . - As described in the 
section ANALYTICAL PROCEDURE the local rates of heat input to the fluid in the feed- 
line were determined from changes in the measured local feedline temperatures with re- 
spect to time. These measured local temperatures, however, were outside metal sur- 
face temperatures. 

Because of the nonsteady nature of these flow tests, the local feedline radial wall- 
temperature gradients were not constant during the experiments. In determining the 
heat-transfer rates, changes in the local radial wall-temperature gradients with time 
were not accounted for. 

The effect of changes in these gradients on the calculated rates of heat transfer was 
investigated. For the tests reported herein, the maximum rate of change of the differ- 
ence between the local measured outside surface temperature and the local average wall 
temperature was estimated to be less than ±3° R per second (5/3 K/sec); that is, 
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For this rate of change, the resulting maximum error in the calculated heat input rate 
was approximately ±5 percent. The maximum relative error in the calculated densities 
was also determined to be about ±5 percent. And the difference between calculated and 
measured densities in figure 7 may be as much as ±5 percent due to possible errors in 
the estimated heat inputs. 
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